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» . Abstract 

c0 , 

Class of spherically symmetric Stephani cosmological models is 
examined in the context of evolution type. It is assumed that the 
t-h ■ equation of state at the symmetry center of the models is barotropic 

(p(t) = ap(t)) and the function k(t) playing role of spatial curva- 
OO | ture is proportional to Stephani version of the Friedmann-Robertson- 

■ Lemaitre- Walker scale factor R(t) (k(t) = j3R(t)). 

Classification of cosmological models is performed depending on 
different values and signs of parameters a and (3. It is shown that 
OQ ' for < (hyperbolic geometry) dust-like (a = 0) cosmological model 

<^ • exhibits accelerated expansion at later stages of evolution. 

J> \ The Hubble and deceleration parameters are defined in the model 

k^J ' and it is shown that the deceleration parameter decreases with the 

j_j ■ distance becoming negative for sufficiently distant galaxies. 

a ■ 



Redshift-magnitude relation m(z) is calculated and discussed in 
the context of Snla observational data. It is noticed that the most 
distant supernovae of type la fit quite well to the relation m(z) calcu- 
lated in the considered model (Hq = 65 km/sMpc, f2o < 0.3) without 
introducing the cosmological constant. 

It is also shown that the age of the universe in the model is longer 
than in the Friedmann model corresponding to the same Hq and f2o 
parameters. 

PACS numbers: 98.62.Py, 98.80.-k, 98.80.Es, 98.80.Hw 
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1 Introduction 



In recent years strong observational evidencies appeared supporting the claim 
that the universe is accelerating its expansion. The observations are based 
on supernovae of type la (Snla) used as standard candles in order to cal- 
culate the distance to them pQ. Since in standard Friedmann-Robertson- 
Lemaitre- Walker (FRLW) cosmological models, with nonrelativistic matter 
and radiation, acceleration is not possible, some people argue that the uni- 
verse must be dominated by cosmological constant or by some other smooth 
exotic component influencing negative pressure. There are many different 
physical scenarios trying to explain existence of such exotic form of matter 

[21 El SI El E]. 

In the present paper we propose different approach, in which exotic mat- 
ter is not necessary to drive accelerated expansion. Instead we relax the 
assumption of homogeneity of space leaving the isotropy with respect to one 
point. In other words we replace the usual Cosmological Principle by Coper- 
nican Cosmological Principle [7] or some generalized version of the Coperni- 
can Principle stating that our Galaxy is at the centre of the Universe. From 
philosophical point of view this assumption seems to be very artificial but 
our intention is to show that even such an exotic model does not contradict 
observations [Sj. 

We shall deal with spherically symmetric Stephani models [21 QUI EH E21 
with the metric given by 

2 
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The functions C(r),k(r) and R{j) are not all independent but are related 
with each other with the help of expressions 
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following from the Einstein equations. ( )' denotes a derivative with respect 
to r. Note that in the considered spherically symmetric Stephani models and 
in the given coordinate system, the energy density p(r) is uniform, while the 
pressure p(r, r) is not and depends on the distance from the symmetry center 
placed at r = 0. This is the reason why in such models barotropic equation of 
state (i.e. of the form p = p(p)) does not exist. If we, however, assume some 
relation between p{r) and p(r, r), this could allow us to eliminate one of the 
unknown functions, e.g. C(t). Hence we are left with two unknown functions 
k{r) and R{r). The first one k(r) plays a role of spatial curvature index, 
while the second one R(t) is Stephani version of FRLW scale factor. As far 
as in FRLW models the curvature index is constant (k = 0, ±1) and does not 
change in the process of evolution of the universe, in Stephani models this is 
not the case. Due to matter transfer, caused by the gradient of pressure, the 
curvature index may change with time [15J. Similarly as in FRLW cosmology, 
where the value of k must be determined from observations, in Stephani 
models the choice concerns the whole function k(r). As regards the scale 
factor R(t), some arbitrariness is connected with the choice of time parameter 
t. We remove this arbitrariness by assuming that the time parameter is 
a proper time of an observer placed at the symmetry center of the model 
(r « 0). 

In the next section we determine spherically symmetric Stephani cosmo- 
logical model by assuming barotropic equation of state in the neighbourhood 
of the symmetry center and by assuming a relation between the scale factor 
R{t) and arbitrary function k(r). 

In Section 3 we perform classification of possible cosmological models 
depending on two parameters: a-parameter appearing in the equation of 
state, therefore determining kind of matter filling up the universe, and (3- 
parameter responsible for the global geometry of the model. 

In Section 4 we focus our attention on one model corresponding to vanish- 
ing a-parameter (dust-like matter) and to negative /3-parameter (hyperbolic 
geometry). The model exhibits accelerated expansion at later stages of evo- 
lution. 

Section 5 is devoted to definition of Hubble and deceleration parameters 
in spherically symmetric Stephani models. We show that in the considered 
model, in spite of the fact that local observations of galaxies give positive de- 
celeration parameter, in the case of more distant ones, negative deceleration 
parameter can be observed. 

In Section 6 we discuss the redshift-magnitude relation in the model and 
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fit our theoretical curves to the supernovae la observational data. 

In Section 7 we compare age of the universe in the considered Stephani 
model to the age in FLRW model determined by the same H , Q and a 
parameters. 

In Section 8 we summarize the results. 



2 Determination of the model 

Now we consider an observer placed at the symmetry center of the spherically 
symmetric Stephani Universe (r ps 0). All our physical assumptions will 
concern his neighbourhood. 

First of all we assume that locally, matter filling up the Universe fulfills 
a barotropic equation of state of the standard form 

p(r ~ 0, r) = cxp(r), (6) 

where a is some constant (a > —1). For example a = corresponds to dust 
and a — 1/3 to relativistic matter. For negative a pressure becomes negative 
and we regard an appropriate matter to be in exotic form. An extreme value 
of a (a — — 1) corresponds to cosmological constant. 

Substituting (4) and (5) into (6) we find an explicit form of the function 
C(r) 

C ( T ) = i2 3(l+a)/2( r )' ( ? ) 

where A is some constant. Choosing a curvature function k{r) in a form 

k(r) = (3R(r) (8) 
determines uniquely the function F(t), which now reads 

F(t) = - , 1 } . (9) 

c ^A 2 - (3R 2 +^{t) 

(3 is a constant, which can be either negative, zero or positive. 

The condition (8) is of course not the most general one. The simplicity 
was the only criterion for choosing it. However, one can easily show that any 
condition of the form k(r) = (3R' 7 (t) (0 < 7 < 2) gives similar results to that 
ones obtained in the present paper. 
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Now we are left with the time parameter r which is not uniquely de- 
termined so far. In FRLW models the time parameter is usually chosen 
globally to be a proper time of any comoving observer. In spherically sym- 
metric Stephani models usually different choice is realized. Since there exists 
only one symmetry center in the models it seems natural to distinguish just 
this point. Hence the time parameter is chosen in a special way to be a 
proper time of an observer placed at the symmetry center. It is not a proper 
time of any other observer. To do that we put [16] 



F{t)IKt)± 



vijyr) 

R{r) 



(10) 



what together with (8) and (9) gives a dynamical equation for the Stephani 
version of the scale factor R(t) 

or more familiar form resembling a Friedmann equation 



fdR(t)\ 2 8ttG 



cdt I + *<*> = < 12 > 
where the energy density p(t) is calculated from (4) 

and C a is some constant. New time parameter is now denoted by t and will 
be regarded as a cosmic time. Solution of the equation (11) could give us a 
full information about evolution of the model. What we are interested in is 
how the type of evolution is influenced by the numerical values of parameters 
of a and (3. Note that Barrett and Clarkson [15J also discussed two parameter 
class of spherically symmetric Stephani models but of different type. Their 
models do not admit dust-like (a = 0) equation of state at the symmetry 
center. 



3 Classification of models 

Qualitative classification of spherically symmetric cosmological models de- 
pending on different values of parameters a and j3 can be performed without 
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solving the equation (11). For this purpose we use a standard method [TTJI 
consisting in treating the evolution equation (11) as the energy conservation 
principle for some dynamical system with kinetic energy (dR/cdt) 2 and po- 
tential energy 



R l+3a 

Then the energy conservation principle reads 

2 

cdt 



^1 +V(R) = Q. (15) 



Now, illustrating the potential energy V(R) on a diagram and keeping in 
mind that (dR/cdt) 2 must be positive we deduce qualitative behaviour of 
the scale factor R(t) for different values of a and (3. We do not consider 
the case (3 = since it corresponds to FRLW Einstein-de Sitter model. All 
other models split naturally into two cases depending on the sign of the (3- 
parameter. Moreover in the frame of each case five subcases can be pointed 
out, corresponding to the following ranges of the a-parameter: 

a) a > —1/3, 

b) a = —1/3, 

c) -2/3 < a < -1/3, 

d) a = -2/3, 

e) -1 < a < -2/3. 

Appropriate behaviours of the potential energy V(R) and corresponding 
evolution of the scale factor R(t) are illustrated in Figs 1 and 2. 

Some of the above presented cosmological models were already discussed 
in detail in literature. Especially much attention was paid to the model 
corresponding to a = —1/3 (Figs lb, 2b) due to its particular simplicity 
(scale factor can be calculated explicitly) [T6], [T9] . In fact, some other models 
can also be expressed in an analytic way. In the present paper, however, we 
shall concentrate on the case, which mathematically is not so simple, but 
from physical point of view is very interesting because it relates to the Snla 
observations suggesting accelerated expansion of the Universe at the present 
epoch pp. The model we want to discuss is presented in Fig. 2a. 
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1 C) p>0, -2/3«x<-l/3 




1 d) P>0, a=-2/3, P<A 2 




Figure 1: Potential energy V(R) of associated dynamical system (left di- 
agrams) and corresponding evolution of the scale factor R(t) (right dia- 
grams) for positive /5-parameter (closed geometry) and different values of 
a-parameter. 



2 a) p<0,a>-l/3 



t 






Figure 2: Potential energy V(R) of associated dynamical system (left di- 
agrams) and corresponding evolution of the scale factor R(t) (right dia- 
grams) for negative /3-parameter (open geometry) and different values of 
a-parameter. 



4 Accelerating Universe 



It is well known that in FRLW cosmology accelerated expansion of the uni- 
verse may be only obtained in the case when the universe is filled with some 
exotic form of matter influencing negative pressure, e.g. cosmological con- 
stant. However, if we look at Fig. 2a we notice that in the spherically 
symmetric Stephani model the universe can exhibit accelerated expansion 
also in the case of vanishing or positive pressure at least at the symmetry 
center. Strictly speaking accelerated expansion takes place for every value 
of a at least at later stages of evolution, provided that (3 < (hyperbolic 
geometry). But it seems to us that the model corresponding to vanishing a 
is the simplest one, moreover it describes locally (in the neighbourhood of 
the observer located at the symmetry center) pressureless matter-dominated 
universe which is in agreement with local observations. 

So putting a = we find the energy density pit) and the pressure p(r, t) 
to be equal to 

, , 3c A 1 , _ „ . 

" (t) = W RHt) (16) 

and 

P ^ t) = ^G Pr ' 2m - (17) 
Note that indeed at the symmetry center (r m 0) pressure vanishes, however, 
due to hyperbolic geometry (5 < it effectively becomes negative at larger 
distances from the symmetry center. So in spite of the fact that the energy 
density scales with the expansion as ordinary nonrelativistic matter p oc 
R~ 3 ), due to negativeness of pressure, for r > 0, it is exotic. Stephani 
version of the Friedmann equation now reads 

(dR\ 2 Qr) A 2 

(18) 

and unfortunately cannot be solved in an elementary way. However in two 
extreme cases corresponding to small and large R exact solutions can be 
easily found. They are 

'3 . \ 2 / 3 



and 



R(t) = (-Ac) t 2/3 for small R, (19) 



R(t) = ^c 2 \P\t 2 for large R. (20) 



9 



In the first case the solution is of Einstein-de Sitter type, while in the second 
case the evolution is evidently inflationary of power law type. It means 
that in our model the universe starts as decelerating and finally ends up 
as accelerating one. In the simplest FRLW cosmological models with one 
component fluid filling up the universe such behaviour is not possible. In 
the present model, inspite of the fact that formally it is of one component 
type (dust-like matter) the curvature term is non-trivial (see the condition 
(8)) and simulates the existence of some exotic fluid driving the power-law 
inflation at later epoch. 

In the next Section we discuss some observational parameters in the 
model. 



5 Observational Parameters 

Observational parameters which describe kinematics of the evolution of the 
universe are Hubble parameter H and deceleration parameter q. In FLRW 
models they are constant on time slices due to homogeneity and isotropy. It 
does not to be the case in spherically symmetric Stephani models, in which 
homogeneity does not take place. 

Following Ellis [20] we define H and q in a general case by 

rr 1^ (dlY 2 (d 2 l\ 1 

where r is here proper time measured along the particle world line, and I is 
some representative length along the particle world-lines. 

In FLRW cosmological models the above expressions go over into the 
standard ones 

R RR 

H=—, q = ^, 22 

R R 2 K J 

where R is a scale factor entering the space part of the line-element 

dr 2 
1 — kr 



ds 2 = c 2 dr 2 - R 2 (t) ( + r W . (23) 



In our spherically symmetric Stephani models the line element is given by 

ds 2 = - (dr 2 + r 2 dn 2 ) . (24) 

V 2 (r,t) V 2 (r,t) \ J K ' 
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Hence neither t is proper time of any comoving observer, nor R(t) is a true 
scale factor (uniquely determining the length I). In order to calculate H and 
q we have to rewrite the line element using global time parametrization 

ds 2 = c 2 dr 2 _ R 2 (r , r) ^ r 2 + r 2 dQ 2^ ^ 

where we defined 

*(r>T) = ^ t{r ; T \ - ( 26 ) 

V 1 K 1 ' V(r,t(r,r)) v ; 

Now the definitions for H and q read 

1 dR(r, r) 1 ^(r, r) 

i2(r,r) rfr ' 9 H 2 R(r,r) dr 2 ' 1 J 

Note that generally and g are not spatially constant anymore, what should 
not be surprising in non-homogeneous universe. It can be shown, however, 
that it is possible to choose time parameter with the aid of which defined 
Hubble parameter does not depend on r. It turns out that the time parameter 
t in the line element (24) is just such a parameter, i.e. 

1 dRjr, r) _ 1 dR(t) 
R{r,r) dr R(t) dt 1 ' 

In other words with this time parameter the expansion is the same at every 
point of time slice of the model. Unfortunately it does not happen in the 
case of q parameter which is evidently coordinate dependent |21j . Explicit 
calculation gives 

q(r, t) = -V(r, t) + V(r, t) - 1 = -V(r, t)q(t) + V(r, t) - 1 (29) 

or 

q(r,t) = q(t) + ±f3r 2 R(t)(l + q(t)), (30) 



where we defined 



»<«> - -mm- (31) 



for deceleration parameter in neighbourhood of the symmetry center. 
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Since (3 < in our model, the deceleration parameter q(r, t) decreases 
with increasing distance to the observed galaxy. Physically it means that in 
our non-homogeneous universe acceleration of the expansion increases with 
distance. 

In order to see how the deceleration parameter q(r,t) depends on the 
redshift we express (3, r and R(t) in terms of observational parameters H , Q 
and z, where Qq is energy density parameter at present epoch and z is a 
redshift of a galaxy placed at some point determined by comoving coordinate 
r. Since in the considered Stephani models energy density pit) is constant 
in space, energy density parameter Q can be defined similarly as in FLRW 
models, i.e. 

m = -^r y (32) 

Pcr{t) 

where 

Per® = j^H\t) (33) 

is critical energy density defined as actual energy density in the model with 
flat geometry (/3 = 0). 

Expressing the comoving coordinate r in terms of observational parame- 
ters is essential for our purposes and can be performed in a standard way by 
putting the infinitesimal space-time interval ds 2 connecting two light events 
equal to zero (we choose coordinate system in such a way that 9 and are 
constant along light ray path from galaxy to observer). Hence we get 

cdt = -R(t)dr (34) 

or equivalently 

-dr (35) 



y/A 2 - (3R 2 +* a 

after getting rid of t variable using (11). Constants (3 and A 2 can be found by 
writing down Friedmann equation (11) or (12) for the present epoch (t = t ) 
and by making use of the definition of Q parameter: 

13 = \r H 2 (Q -1), (36) 



c 



A 2 = \R 3+3a H 2 Q . (37) 



c 2 ' 
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Substituting these expressions into (35), next integrating both sides and 
changing integration variable we get integral formula for r 



„(3a— 1)/2j„ 

(38) 



RoH 0r J Rq ^ + (i _ n )x^ 

R is a scale factor at the epoch when the light ray was emitted by the observed 
galaxy. Unfortunately integration cannot be performed explicitly either in 
general case or even in the special case corresponding to a = 0. We compute 
it approximately assuming that Qq is small with respect to the ratio R/Rq. 
Physically it means that in order to keep the approximation valid, the larger 
distances r are calculated, the smaller energy density parameter fl has to 
be taken into account. Then we skip the first term under the square root in 
(38) and the integration can be carried out explicitly giving 



2c 



RoH y/l — Q 
Together with the relation [16] 

R z+1 




(39) 



R V{r,t) 



(40) 



we obtain system of algebraic equations which can be solved with respect to 
R and r 

* = ™ (41) 



(z + 2) 2 ' 

cz 



RqHq^I — VLq 



(42) 



Substituting (28a) and (33) into (24b) we find present value (t = t ) of the 
deceleration parameter of the region of the universe where a galaxy with the 
observed redshift z is placed 

q(z,t ) =q -^{l + q )z 2 . (43) 
q Q is a local value of the deceleration parameter. 
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In this model, in spite of the fact that local observations of galaxies give 
positive value of q , in the case of more distant ones (z > y^o70~+~9o)) 
negative deceleration parameter is measured. 

More detailed description of the model from the point of view of observa- 
tions requires derivation of coordinate independent relation between obser- 
vational parameters. Redshift-magnitude relation is just such a relation and 
finding it will be our next task. 



6 Redshift-magnitude relation 

In cosmological models which are not homogeneous and isotropic one of the 
methods allowing for finding redshift-magnitude relation is based on formal- 
ism proposed by Kristian and Sachs [22] and Ellis and MacCallum [23] leading 
to relation in form of power series around the observer's position. But this 
method is not very effective if large redshifts are taken into account. 

Due to spherical symmetry of the model it is possible, however, to find 
closed formula for apparent bolometric magnitude m of galaxy, as a function 
of redshift z in the case of small ^-parameter. The problem consists then 
in explicit calculation of the luminosity distance Dq. 

The formula for the luminosity distance in our spherically symmetric 
model is the same as in FLRW models and reads 



D = R r(z + l). (44) 
Making use of earlier derived expressions for /3, r and R we find 

For small redshifts (z < 1) the expression goes over into usual distance- 
redshift relation known from Friedmann cosmology (we remind that the above 
formula is valid only for small values of Qq) 

Do « ^-z. (46) 

However, for larger redshifts it grows faster than in FLRW models, where 

D ° = IFT^^oz + (fi - 2)( V /O z + 1 - 1)], (47) 
_H "0 
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and is in favor of Snla observations. 

Finding general formula for D requires numerical integration of (38). 
Before doing that we expand Dq (z) in our model into Taylor series with 
respect to the z-variable. Up to the second order we get (for arbitrary fl 
and a) 

D (z) = ±. |z + j [4 - 3n (a + . (48) 

Comparing with the analogous expansion of the Mattig formula [25], which 
up to the same order in the dust filled universe (a = 0) reads D (z) = 
(c/H )[z + (2 — Qo)z 2 /4], we notice that for arbitrary Q redshift-magnitude 
formula in the Stephani model indeed grows faster than in the Friedmann 
model at least for relatively small values of z, for which the Taylor expansion 
is valid. 

Since supernovae of type la favoring accelerated evolution of the Universe 
were discovered at relatively large redshifts z > 0.5 full formula for the 
luminosity distance is required valid not only for small z. Hence Taylor 
expansion calculated up to low order is not satisfactory and we must calculate 
Dq(z) numerically. 

The results of this numerical calculation are presented in figures 3 and 
4. First (figure 3), we compare the redshift magnitude relation obtained in 
both Friedmann and Stephani spherically symmetric models determined by 
the same cosmological parameters H , Q and the a-parameter responsible for 
the form of matter filling locally the universe. We realize that for relatively 
large redshifts the curves m(z) in both models diverge and the relation m(z) 
in the Stephani model looks like the appropriate relation in the Friedmann 
model but with the cosmological constant. In other words nonhomogeneity 
in the spherical symmetric Stephani model in some sense simulates existence 
of the A-term. 

In Fig. 4 the redshift-magnitude relation in the Stephani model is com- 
pared with some observational data from the Supernova Cosmology Project 
[26] (see also [27]). We notice that the most distant supernovae of type la 
fit quite well to the dust-like (a = 0) Stephani cosmological model with the 
energy density parameter Q < 0.3. The values of the Hubble constant H 
and the absolute magnitude M are the same as in [26l [28] . 

In the next section we compare numerically the age of the universe in the 
considered model and in FLRW cosmological models corresponding to the 
same value of H and fl parameters. 
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Figure 3: Redshift-magnitude relation m(z) in Friedmann cosmological 
model (lower line) and in the spherically-symmetric Stephani model (upper 
line). In both models the same cosmological parameters are chosen: Hq = 65 
km/sMpc, M = -19.5, fl = 0.3. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Oo 

Figure 5: Age of the universe in spherically symmetric Stephani cosmo logical 
model (upper line) and in FLRW cosmological model (lower line) as a func- 
tion of energy density parameter Q . It is assumed that a = and H = lOO/i 
km/sMpc, where 1/2 < h < 1. 



7 Age of the universe 

We find the age of the universe in the considered models by direct integration 
of Stephani version of the Friedmann equation (11) 

1 ^° R (l+3a)/2 dR l 1 X {^)/2 dx 



Analogously calculated age of the universe in FLRW models reads 

p 1 } x^ 1+3a ^ 2 dx , x 

to F = -FT , — • 50 



We easily realize that the age in Stephani models is longer than in FLRW 
ones corresponding to the same values of the parameters Hq, Qq and a. For 
a = and Hq = 100/i km/sMpc (1/2 < h < 1) and on the same diagram 
(Fig. 5) we plot the age of the universe, to S and to F , in both models as a 
function of energy density parameter Qq- 
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Figure 6: Age of the universe in spherically symmetric Stephani cosmo logical 
model as a function of a and Qq. 



In Fig. 6 we plot the age to (a,Qo) of the Universe in Stephani model 
as a function of two parameters: a - determining equation of state at the 
symmetry center, and Qq - energy density parameter. We notice that the 
age increases with decreasing both a and Qq (provided that a > — 2/3). In 
Fig. 7 we present sections of the plot t s (a, Qq) corresponding to the values 
of a =-2/3, -1/3, 0, 1/3, 2/3, 1. 

8 Conclusions 

In the paper we considered a class of spherically symmetric Stephani cosmo- 
logical models parametrized by two parameters a and (3. The a-parameter 
determines the form of barotropic equation of state in the neighbourhood of 
the symmetry center of the models and /3-parameter is related to the cur- 
vature index. We performed qualitative classification of evolution types of 
models depending on different values of a and (3. We showed that for nega- 
tive P the universe exhibits accelerated expansion independently of a, what 
does not take place in FLRW models. We focused our attention on the model 
corresponding to a = 0. In this model pressure is negligible in the neighbour- 
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Figure 7: Age of the universe in spherically symmetric Stephani cosmo- 
logical model as a function of f2 for different values of a = —2/3, — 
1/3, 0, 1/3, 2/3, 1. 



hood of the symmetry center (the so-called dust-like model) and the model 
looks like of Einstein-de Sitter type at the early stage of evolution. At later 
epoch the evolution is inflationary of power law form. This result seems to be 
interesting because we showed that existence of "exotic" matter (influencing 
negative pressure) is not necessary (at least locally) to drive accelerated ex- 
pansion of more distant regions. We also showed that acceleration becomes 
larger while increasing the distance to the galaxy. 

We derived redshift-magnitude formula in the model and showed that 
for large relatively large redshifts (z ~ 1) apparent bolometric magnitude 
m(z) grows faster than in corresponding FLRW model what is in favor of 
Snla observational data. We found numerically that a good fit is obtained 
for H = 65 km/sMpc and Qo < 0.3. At this point it should be noted that 
the good fit of m(z) to the SNIa observational data in spherically-symmetric 
Stephani cosmological model has been already obtained by Dabrowski and 
Hendry [28J. However in their model they assumed a string-like (a = —1/3) 
equation of state at the center of symmetry. 

Finally we found that the age of the universe in the considered model is 
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remarkably larger than in isotropic and homogeneous models. 

Concluding we would like to stress that we do not claim that our sur- 
rounding universe is of Stephani type. We wanted to show that there exist 
different from isotropic and homogeneous cosmological models which give 
the same observational evidencies as FRWL models with the cosmological 
constant. In these models the deviation from homogeneity could be so small 
that locally (up to several hundreds of megaparsecs)it would be undetectable. 
However its existence could show up at larger distances. 

In the paper we considered only one fluid component filling up the uni- 
verse. It is of course oversimplification and is hardly acceptable especially in 
the context of negativness of pressure at larger distances. It can be shown, 
however, that adding usual dust component into the model do not change 
relations (16, 17) significantly (only a constant in (16) will be different) and 
the qualitative result will be the same. 

It seems that these two different fluid components should effect in the 
form of the power spectrum of the CMB anisotropies. Calculation of this 
power spectrum in the proposed model will be our next task. 

We thank M.P. Dabrowski for valuable discussions. 
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